Thermally stable carbon-related centers in 6H-SiC: 
photoluminescence spectra and microscopic models 
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Recent ab initio calculations [Mattausch et al, Phys. Rev. B 70, 235211 (2004)] of carbon 
clusters in SiC reveal a possible connection between the tricarbon antisite (C3)si and the U pho- 
toluminescence center in 6H-SiC [Evans et at, Phys. Rev. B 66, 35204 (2002)]. Yet, some of the 
predicted vibrational modes were not observed experimentally. Here we report experiments which 
indeed confirm the existence of a low-energy mode for the U-center (as well as for the HT3- and 
HT4-centers with spectral details similar to the U-center). We calculated the isotope splitting for 
the (C3)si-defect and found near-perfect agreement with our data. In addition, we discuss the car- 
bon di-interstitial (C2)Hex as a model for the Z- and HT5-centers. The isotope splitting is also well 
reproduced, but the absolute values of the local mode energies show a discrepancy of about lOmeV. 



PACS numbers: 61.80.-x, 63.20.Pw, 78.55.-m 

Silicon carbide (SiC) is a wide band-gap semiconduc- 
tor especially suitable for high-power, high-temperature 
and high-current applications. Due to its extreme hard- 
ness and low diffusion coefficients, the method of choice 
for the creation of the doping profiles is ion implanta- 
tion. However, implantation causes a severe damage of 
the crystal and a subsequent annealing at high tempera- 
tures is necessary to eliminate the defects — mostly vacan- 
cies and interstitials. Unfortunately, not all interstitials 
and vacancies recombine. They can also cluster and form 
thermally stable aggregates. One example is the well- 
known Dn-centerji which is most likely an aggregation 
of carbon atoms around carbon antisites^*^ 

Recent photoluminescence (PL) measurements^ of 
electron-irradiated 6H-SiC have shown that besides the 
Dn-center other intrinsic defect centers exist which also 
possess localized vibrational modes (LVM's) above the 
SiC bulk phonon energy. These are the P— T cen- 
ters, which were tentatively assigned to carbon split- 
interstitials or di-carbon antisites,^''* and the U-center, 
for which a single LVM at nearly 250 meV was observed. 
The three- fold isotope splitting in ^'^C-enriched samples 
signifies a carbon-related defect, although carbon-based 
defects with such a high vibrational frequency have not 
been known in SiC. Yet, ab initio calculations predict^ 
that such defects may exist. Namely the tri-carbon an- 
tisite (C3)si has several LVM's up to 250 meV. The en- 
ergetically highest mode even reproduces the three-fold 
isotope splitting of the U-center. However, the missing 
phonon replicas in available PL spectra could not be ex- 
plained in Ref.0. Yet, the clear identification of observed 
carbon-related defect centers is important since it may 
substantiate the carbon aggregation phenomenon in SiC 
predicted by earlier theoretical work i^i^i^ Besides, at ele- 



vated temperatures these aggregates can re-emit carbon- 
interstitials^ and thereby facilitate the formation of fur- 
ther defect centers like Dj or the alphabet-lines J 

In this paper we present photoluminescence results for 
^^C-enriched 6H-SiC samples and compare them with ab 
initio calculations. The data show an additional phonon 
replica of the U-center at 150 meV, which is in agree- 
ment with the calculated LVM of (C3)si. Further cen- 
ters, which are similar to the U-center and labeled HT3 
and HT4, were also observed. With this new experimen- 
tal data in hand, the centers are well explained by the 
tri-carbon antisite (C3)si model. The centers Z and HT5 
are also discussed. These centers possess two vibrational 
modes above the SiC bulk phonon spectrum at about 
170meV and 200 meV and also exhibit a three- fold iso- 
tope splitting of the energetically highest mode. Their 
relationship with the carbon di-interstitials is discussed. 
We show that a three-fold isotope splitting with the fre- 
quency ratios of a C— C dumbbell (which is the square 
root of the mass ratios), does not necessarily originate 
from a dumbbell-shaped defect. Both the tri-carbon an- 
tisite and the carbon di-interstitial represent complex de- 
fect configurations. 

In the following, we first present the experimental de- 
tails and results and then turn to the theoretical calcu- 
lations. Theoretical results are then compared with the 
experimental data. 

The investigated samples were n- and p-type 6H-SiC 
with typical doping levels in the range 10^^ — 10^® cm^'^. 
They included both epitaxially grown layers and Lely 
grown crystals. The p-type specimens showed strong 
donor-acceptor pair-luminescence. To produce the de- 
fect centers the samples were electron irradiated with 
electron energies below and above the silicon displace- 



2 



TABLE I: ZPL's and LVM's of the centers in 6H-SiC. 
center ZPL (nm) ZPL (eV) LVMl (meV) LVM2 (meV) 



z 


512.5 


2.417 


171.6 


202.0 


HT5 


509.8 


2.431 


169.0 


201.7 


HT3 


492.9 


2.515 


151.0 


246.3 


U 


525.0 


2.361 


151.6 


246.6 


HT4 


546.4 


2.268 


149.9 


244.7 



ment threshold at about 250 keV. The irradiation was 
performed using an ion-free transmission electron micro- 
scope. To allow for an identification of the defects via the 
isotope splitting of the phonon replicas, samples with a 
^^C enrichment of 30% were used in addition to the sam- 
ples with natural isotope abundance. For the PL mea- 
surements the samples were cooled to 7K using Oxford 
Instruments microstats on Renishaw microRaman spec- 
trometers. The centers were excited using laser wave- 
lengths of 325 nm and 488 nm. 

The reported centers were produced by high- 
temperature anneals in the range 1000— 1500°C. The U- 
and the Z-center came up at the lower end of this range 
and persisted in some cases up to 1500°C. The centers 
HT3 and HT4 were observed after 1200°C and 1300°C 
anneals, respectively. While HT3 was always present, 
HT4 appeared only in some samples. The HT5-center 
was observed after the highest temperature of 1500°C. 

All of the optical transitions were excited by both the 
325 nm and the 488 nm wavelength. A major effort was 
put into the deconvolution of the observed spectra, since 
considerable overlap occurred between the spectral de- 
tails. Employing three different techniques we were able 
to unambiguously assign the spectral details to a specific 
zero phonon line (ZPL). One method was by spatial sep- 
aration. Spectra were obtained at a series of points along 
lines through the irradiated region, and the relative in- 
tensities of the ZPL's changed with position. A second 
approach was the comparison of the results after different 
annealing temperatures, which also affected the relative 
intensities of the spectra. The third method was the in- 
vestigation of a wide range of specimens irradiated with 
different electron doses and energies. Further details on 
this approach can be found in Refs.0 andH. 

Figure nji shows the PL spectra of the U- and the Z- 
centers in ^^C-enriched 6H-SiC. Both centers were origi- 
nally discussed in Ref. |^ The three- fold isotope splitting 
of the high-energy mode at 246.6 meV is clearly visible. 
In addition to the single originally observed mode^ a fur- 
ther mode with a lower energy was observed at 151.6 meV 
in samples with natural C abundance (cf . Table P) . In 
isotope-enriched SiC this mode splits in a complex fash- 
ion (cf. Fig.^). The centers HT3 and HT4 possess prop- 
erties similar to those of the U-center. Figure^ shows 
the HT3-center alongside the Z-center in ^^C-enriched 
6H-SiC. Each center also has a high-energy mode around 
245 meV (cf. Tabled exhibiting a three-fold isotope split- 



ting (cf. inset of Fig. QJj for HT3) and a lower-energy 
mode around 150 meV. Some aspects of the HT3- and 
HT4-centers deserve particular mention. First, in the 
sample with natural C-abundance the high-energy mode 
coincides with the ZPL of HT4 at 546 nm (cf. Fig. Ut). 
Clear evidence of a LVM at this wavelength could be de- 
duced from a ^^C-enriched sample which lacks the HT4- 
luminescence and shows a clear three-fold isotope split- 
ting. The second aspect is the small separation between 
the lower-energy mode of HT3 and the ZPL of the U- 
center (cf. Fig. QJ;). Although the LVM could be clearly 
identified in the sample with natural C-abundance, the 
isotope splitting could not be deduced. 

The Z-center is clearly present in all shown spectra. Its 
full details are visible in Fig. It possesses 4 LVM's, 
yet only two of them (Zq and Zd in Fig. ^p) are above 
the SiC bulk phonon spectrum at 172 meV and 202 meV. 
The HT5-center possesses similar properties (cf. TablePl. 
In ^^C-enriched samples the high-energy mode shows a 
three-fold isotope splitting, while the lower-energy mode 
exhibits a complex broadening (cf. Fig. ^). 

The high vibrational energy and the three-fold split- 
ting of the highest mode upon carbon isotope substitu- 
tion suggest that the observed defect centers are carbon- 
related. Hence we focused our theoretical modeling on 
carbon-interstitial and carbon-antisite complexes. We 
employed density-functional theory with the local den- 
sity approximation (LDA) as implemented in the soft- 
ware package FIII96SPIN.^ Supercells with 216 sites for 
3C-SiC and 128 sites for 4II-SiC and a plane wave basis 
set with a cut-off energy of 30 Ry were used. The Bril- 
louin zone was sampled by a single F-point for the large 
3C-SiC supercells. For 4H-SiC a 2 x 2 x 2 Monkhorst- 
Pack mesbiS was used except for the calculation of the 
LVM's, where we used only the F-point. Due to the 
large supercells the defect molecule approximation was 
employed for the analysis of the LVM's. In this approach 
the calculation of the dynamical matrix is restricted to 
the defect and its nearest-neighbor atoms (cf. Ref. 3 for 
details). We verified that this approximation affects only 
the lowest-energy LVM's close to the bulk phonon bandit 

The candidates for the centers described above must 
possess two properties: they must exhibit a three-fold iso- 
tope splitting of the highest mode in ^^C-enriched mate- 
rial and they must be thermally stable. Two possibilities 
are the carbon di-interstitial (C2)hcx and its structurally 
similar configurations in 4II-SiC and the tri-carbon anti- 
site (C3)si. The calculations were performed for 3C- and 
4H-SiC. A close agreement of the results for the two poly- 
types ensures the transferability to the more complex 611- 
SiC. The defects are depicted in Fig. for 3C-SiC. The 
thermal stability is guaranteed by high dissociation ener- 
gies, which are 4.8 eV (5.8 eV) for (C3)si in 3C-SiC (4H- 
SiC) and 4.8 eV for (C2)Hex in 3C-SiC. In 4H-SiC, various 
configurations exist for the di-interstitial with dissocia- 
tion energies between 5.1 eV and 5.5 eV. The most stable 
defect is (Csp)2,hh with two carbon split-interstitials at 
neighbouring sites. The interstitial carbon atoms relax 
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FIG. 1: Spectra of 6H-SiC samples (taken at T = 7K with a laser wavelength of 325 nm) with a ^^C-enrichment of 30% and 
natural isotope abundance showing the U-, Z-, HT3- and HT4-centers. The shaded regions indicate the bulk phonons belonging 
to the marked ZPL. (a) The U- and the Z-center in a ^^C-enriched sample. The inset displays the broadened lower-energy mode 
of the U-center and the ^^C— ^^C-peak of the energetically highest Z-mode. (b) The HT3- and the Z-center in a ^^C-enriched 
sample. All modes of the Z-center are visible. The inset displays the isotope splitting of the modes Zd and Zc as well as the 
three-fold splitting of the highest HT3-mode. (c) The centers Z, U, HT3 and HT4 in a sample with natural isotope abundance. 
The overlap of the ZPL of HT4 and the high-energy mode of HT3 is clearly visible, as well as the narrow separation between 
the lower-energy mode of HT3 and the ZPL of the U-center. 




FIG. 2: The tri-carbon antisite (C3)si (left) and the carbon 
di-interstitial (C2)hcx (right) in 3C-SiC. 



towards each other, so that a configuration similar to 
(C2)hcx in 3C-SiC forms. All these defects are neutral 
for practically all positions of the Fermi-level. 

Both defects shovif a distinct vibrational pattern. The 
frequencies of the LVM's in 3C- and 4H-SiC are listed 
in Table ^ Since the defects possess mirror symmetry 
(which is in 4H-SiC limited to the defect molecule), the 
modes are marked as symmetric (g) or antisymmetric (u) 
vibrations. This is important for PL experiments, since 
the bound exciton mostly interacts with a symmetric but 
not with an antisymmetric vibration. It is clearly visible 
from Tablelnlthat the LVM's of the defects are practically 
independent of the polytype and the involved cubic or 
hexagonal sites. We therefore expect that our results for 
the 3C and 4H polytypes are transferrable to 6H-SiC. 

The di-interstitial possesses four LVM's above the SiC 
bulk phonon spectrum. The LVM's at 160 mev and 
190 meV are symmetric vibrations (cf. Table llT|l . The 
highest mode shows a clear three-fold isotope splitting, 
while the lower-energy modes show a complex broaden- 



ing or splitting (cf . Fig. . The relative height of the 
peaks within a mode results from the various possibili- 
ties of substituting a defect's atom with ^'^C assuming an 
enrichment of 30%. The peak heights of different modes 
are independent and cannot be compared. Considering 
the two symmetric modes of (C2)hcx, the spectrum is 
similar to that of the Z- and the HT5-center. Besides, as 
in 4H-SiC various different defect configurations involv- 
ing the two cubic and the hexagonal sites are possible 
in 6H-SiC, allowing for different exciton binding energies 
(and therefore different ZPL's) with comparable phonon 
replicas. However, the calculated values are systemati- 
cally about lOmeV lower than the measured frequencies. 
This finding is especially noteworthy since the LDA tends 
to overestimate the calculated frequencies. 

The tricarbon antisite (C3)si has five LVM's above the 
bulk phonon spectrum. Symmetric modes are the mode 5 
at 250 meV and the mode 3 at 150 meV. All other modes 
are antisymmetric (cf. Table . The isotope splitting 
of the defect in 3C- and 4H-SiC is displayed in Fig. 
alongside the measured frequencies for the U-center (the 
centers HT3 and HT4 have similar phonon replicas, cf. 
Table P) . The mode 5 shows a clear three- fold isotope 
splitting as one would expect for a simple carbon dumb- 
bell vibration. This simple pattern results from the lo- 
calization of the vibration onto the atoms Ci and C2 of 
the triangular defect, whereas the neighboring atoms are 
practically not involved. In contrast, the neighbor atoms 
do participate in the vibration of the mode 3. Conse- 
quently, this mode exhibits a significant broadening due 
to the ^"^C-enrichment. The measured values of the U-, 
HT3-, and IIT4-centers are in excellent agreement with 
the calculated absolute values of the symmetric modes of 
(C3)si. The isotope splitting of the centers is also well 
reproduced. The antisymmetric modes should be sup- 
pressed in photoluminescence as discussed above. In 611- 
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TABLE II: Calculated LVM's in meV of the carbon di-interstitials (C2)hcx and the tri-carbon antisite (Cajsi- In 4H-SiC, 
further di-interstitials exist with similar geometrical configuration and similar vibrational properties (see Ref. Q)- The letters g 



and u 


denote symmetric and antisymmetric vibrations 


respectively. 








LVM 


(C2)hcx (3C) 


(Csp)2,hk,cub (4H) 


(ap)2,hh (4H) 


(C3)si (3C) 


(C3)si,k (4H) 


(C3)si,h (4H) 


1 


132.5 (u) 


128.0 (u) 


134.9 (u) 


118.6 (u) 


119.0 (u) 


119.0 (u) 


2 


160.4 (g) 


161.3 (g) 


159.9 (g) 


129.8 (u) 


130.2 (u) 


130.5 (u) 


3 


167.3 (u) 


167.1 (u) 


168.1 (u) 


148.7 (g) 


154.0 (g) 


154.2 (g) 


4 


184.3 (g) 


189.1 (g) 


191.5 (g) 


181.3 (u) 


182.3 (u) 


182.6 (u) 


5 








248.5 (g) 


254.9 (g) 


255.7 (g) 



(C,)Hex [3C-SiC] 



(Csp)2,hh [4H-SiC] 




iJ 



1 



(b) 
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200 100 
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Frequency of LVMs (meV) 
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FIG. 3: (Color online) Calculated LVM's of (a) the carbon di-interstitials and (b) the tri-carbon antisite in 3C and 4H-SiC 
with a ^^C-enrichment of 30%. The black peaks indicate symmetric modes, the red (gray) peaks are antisymmetric vibrations. 
The phonon replicas of the Z-center (a) and the U-center (b) including isotope shifts are indicated by the blue (shaded) lines. 



SiC the defect can be located at the cubic sites fci and 
/c2 or at the hexagonal site h, providing three different 
defect configurations. As observed for 4H-SiC, the dif- 
ference of the vibrational patterns for the different sites 
shall be negligible. Thus we believe that (C3)si should be 
a good model for the centers U, HT3 and HT4 in 6H-SiC. 

In conclusion, we presented thermally stable photolu- 
minescence centers in 6H-SiC with localized vibrational 
modes up to 250 meV. All centers show a clear three- 
fold isotope splitting of the highest mode with the fre- 
quency ratio typical for a carbon-carbon dumbbell vi- 



bration. Based on the vibrational spectra obtained from 
ab initio calculations, the di-interstitial (C2)hcx is con- 
sidered as a model for the centers Z and HT5, although 
a systematic frequency shift of 10 meV is observed. The 
results obtained for the tri-carbon antisite (03)31 are in 
excellent agreement with the U-, HT3-, and HT4-centers. 
This identification demonstrates the possibility of carbon 
aggregates predicted earlier theoreticallyi2i2i4 
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gemeinschaft within the SiC Research Group and the UK 
Engineering and Physical Science Research Council. 
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